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The search for resonant absorption of the Primakoff solar axions by ^^^Tm nuclei have been performed. Such an 
absorption should lead to the excitation of low-lying nuclear energy level: A+^^^Tvci ^^^^Tm* ^^^^Tm +7 (8.41 keV). 
The Si(Li) detector and ^^^Tm target placed inside the low-background setup were used for that purpose. As a result, 
a new restriction on the axion-photon coupling and axion mass was obtained: ^^^(GeV^^) • mA{eV) < 1.36 • 10~^ (90% 
C.I.). In model of hadronic axion this restriction corresponds to the upper limit on axion mass - tua < 191 eV for 90% 
c.l. 
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1. Introduction 

The appearance of an axion in theory is connected with 
the problem of CP-violation in strong interactions. In or- 
"der to solve this puzzle Peccei and Quinn [l[ proposed 
.the concept of the new chiral symmetry U{1). The spon- 
taneous breaking of this symmetry at the energy /a al- 
lows one to compensate CP-violating term of the QCD La- 
■grangian completely. Weinberg 0] and Wilczek Q showed 
that the introduced model should lead to the occurrence of 
"a new pseudoscalar particle. The axion mass rriA appears 
,to be inversely proportional to the Ja value, as well as the 
feffective axion coupling constants with photons (3^7), lep- 
,tons {gAe), and hadrons (gAN)- The model of "standard" 
or PQWW-axion, where the value of Ja was fixed at the 
electroweak scale {/a ~ v^Gi?)"^/^ has been excluded by 
the series of experiments with radioactive sources, reactors 
and accelerators. 

' Two types of the "invisible" axion models retained the 
axion in the form required for the solution of CP-violation 
problem, while suppressing its interaction with matter, 
since the scale of symmetry breaking fA appears to be arbi- 
trary in these models and can be extended up to the Planck 
mass mp w lO^^GeT^. These are models of "hadronic" or 
KSVZ axion [1,0] and the GUT or DFSZ axion The 
axion mass ttia (in eV units) in both models is given in 
terms of tt" -properties: 



ruA = 
where /^r 



1/2 



6.0- 10*^ 



fA 'il + z + w)il + z)' fAiGeV) 
': 93 MeV is the pion decay constant; z 



(1) 



"Tnu/Tiid = 0.56 and w ~ mu/ma = 0.029 are quark-mass 
ratio^j. 

The restrictions on the axion mass appear as a result 
of the restrictions on the coupling constants gA^^ gAe 
and gAN, which are significantly model dependent. The 
hadronic axion does not interact with leptons and ordi- 
nary quarks at the tree level, which results in strong sup- 
pression of gAe coupling through radiatively induced cou- 
pling. Moreover, in some models axion-photon coupling 
may significantly differ from the original DFSZ or KVSZ 
gA~j couplings . The coupling constant gA-y for the "invis- 
ible" axion models is equal to: 
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where a ~ 1/137 is a fine structure constant, E/N is the 
ratio of the electromagnetic and colour anomalies, a model 
dependent parameter is of the order of unity. E/N = 8/3 
for the DFSZ-axion (Ga77=0.74) and E/N = (GA-y7=- 
1.92) in the original KSVZ-axion model. The value of the 
second term inside brackets is 1.95±0.08 and axion-photon 
coupling may be reduced by a factor less than 10^^ in 
axion models in which E/N is close to 2 Q. 

If axions do exist, then the Sun should be an intense 
source of these particles. Axions can be efficiently pro- 
duced in the Sun by the Primakoff conversion of photons 
in the electric field of the plasma. The resulting axion flux 
depends on g\^ and can be detected by inverse Primakoff 
conversion of axions to photons in laboratory magnetic 
fields - [l^l or by the coherent conversion to photons in 



Email addresses: derblnOpnpi . spb . ru (A.V. Derbin) 



^We follow the generally accepted units: GeV for /a and qaj and 
eV for rua- 
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the crystal detectors [l^- 22 1. The expected count rate of 
photons depends on the axion-photon couphng as g\^. 

In this letter we present the results of the search for solar 
axions using the other reaction - the resonant absorption 
by nuclear target. The 7-rays and conversion electrons 
produced by the de-excitation of the nuclear level can be 
registered. The detection probability of the axions is de- 
termined by the product g\^-g\fq which is preferable for 
small QA-y values. 

There are two other possible mechanisms of axion pro- 
duction in the Sun: the reactions of solar cycle and the 
excitation of the low-lying energy levels of some nuclei by 
the high solar temperature. The attempts to detect the 
resonant absorption of quasi-monochromatic axions emit- 
ted in nuclear magnetic transitions were performed in psj - 
3ll |. Astrophysical and cosmological data provide restric- 
tive constraints on axion mass so that acceptable val- 
ues belong to the 10"^ ~ 10~2 eV region [H, i^. The 
results of laboratory searches for the axion as well as the 
astrophysical axion bounds one can find in [33|, l3J| 

2. The rate of solar axions absorption by ^^^Tm 
nucleus 
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The energy spectrum of solar axions produced by Pri- 
makoff effect is parameterized by the following expression 
11, 



35|,llZ|: 



cLEa 



{gA-yf ■ 3.82 • 10 



30 



{Ea) 



exp(£;A/1.103) - 1' 



(3) 



where the value of flux is given in (cm^^s~^keV^^) units, 
value of Ea energy is given in keV units and the value of 
gA-y is given in GeV~^ units. The axion spectrum calcu- 
lated in assumption that gA^ — 10^^° GeV~^ is given in 
Figdl The average energy of axions is 4 keV, but their 
flux becomes negligibly small at the energies above 15 keV. 

As a pseudoscalar particle, the axion should be subject 
to resonant absorption and emission in the nuclear tran- 
sitions of a magnetic type. For our experiment we have 
chosen the ^^^Tm nucleus as a target. The energy of the 
first nuclear level (3/2+) is equal to 8.41 keV, the axion 
flux at this energy is only 7 times less than the maxi- 
mum level (figd]). The 8.41 keV nuclear level discharges 
through Ml-type transition with i?2-transition admixture 
value of 5^=0.11%. Considering the electron conversion 
ratio e/7 = 263 [11] one obtains the relative probability of 
7-ray emission 77 = 1/(1 -|- e/7) = 3.79 • 10"'^. 

The cross-section for the resonant absorption of the ax- 
ions with the energy Ea is given by the expression that is 
similar to the one for 7-ray resonant absorption, but the 
ratio of the nuclear transition probability with the emis- 
sion of an axion {oja) to the probability of magnetic type 
transition {lo^) have to be taken into account. The rate of 
solar axion absorption by ^^^Tm nucleus will be 



Ra 



dtjA 



(4) 



Figure 1: Energy spectrum of the axions produced by Primakoff 
eflfect at the Sun (for 3^7 = 10" 1° GeV). The level scheme of ^^STm 
nucleus is shown in the inset. 



where ctq^ is a maximum cross-section of 7-ray absorption. 
The experimentally derived value of cro-y for ^®^Tm nucleus 
is equal to 0-07 = 2.56 • 10"^^ cm^. A hfetimc of the ^®^Tm 
first excited level is r = 5.89 ns, thus the width of the 
energy level T = 1.13 • 10""^ eV. 

The ujaI'-o^ ratio calculated in the long- wave approxi- 
mation, has the following view [S^, : 



1 



1 





W7 27ra 1 -|- (5^ 



{fio - 0.5)/? -t- - 77 



^ ) .(5) 

P7. 



Here, and pA are the photon and axion momenta, re- 
spectively, = + /i„ ~ 0.88 and ^jl^ = ^ip — ^JLn ~ 4.71 
are isoscalar and isovector nuclear magnetic momenta, [3 
and 77 are parameters depending on the particular nuclear 
matrix elements [ssl, [39| . In case of the ^^^Tm nucleus, 
which has the odd number of nucleons and an unpaired 
proton, in the one-particle approximation the values of (3 
and 77 can be estimated as /? « 1.0 and 77 « 0.5. For the 
given parameters the branching ratio can be rewritten as: 
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In the KSVZ axion model the dimensionless isoscalar 
and isovector couplingconstants g\j^ and g\j^ are related 
to fA by expressions : 



QjA 1 + z + w 



(7) 



2 



and 

where niN ~ 939 MeV is the nucleon mass. The axial- 
coupUng parameters F and D are obtained from hyperon 
semi-leptonic decays with a high degree of precision: F = 
0.462 ± 0.011, D= 0.808 ± 0.006 |4l|. The parameter 
S characterizing the flavor singlet couphng still remains a 
poorly constrained one. Its value varies from S = 0.68 in 
the naive quark model down to S' = —0.09 which is given 
on the basis of the EMC collaboration measurements . 
The more stringent boundaries (0.37 < S < 0.53) and 
(0.15 < S < 0.5) were found in [i^ and accordingly. 
As a result the value of the sum {g an '^9 an) is determined 
within a factor of two, but the ratio ujaI^-^ does not vanish 
for any value of parameter S, in the case of the transition 
caused predominantly by the proton (/3 >0). We will use 
S = 0.5 as reference when we calculate an axion flux for 
KSVZ axion model. 

The values of and g\j^ for the DFSZ axion depend 
on an additional unknown parameter cos^ /? which is de- 
fined by the ratio of the Higgs vacuum expectation values. 
In case of ^^^Tm Ml-transition the value of (w^/w^)^^'^^ 
ratio hes within the interval ~ (0.1-=-2.0)(tt'A/w7)^'^^^- 
The lower and upper bounds of this interval are defined 
by values cos /3 = 0, 1 respectively. 

In accordance with ([3]) and ^ , the rate of axion absorp- 
tion by ^^^Tm nucleus ^ dependent only on the coupling 
constants is (the model- independent view): 

Ra = 104 • gX{g% + g\^)\vAlV'?i' ■ (9) 

Using the relations between g^jy 9\n ^'^'^ axion mass 
given by KVSZ model ([7]), the absorption rate can be pre- 
sented as a function of gA-i and axion mass m^: 

Ra = 4.80 • \Q-^^g\m\{VAlv^f- (10) 

At last, one can use the connection between gA^ and 
given by expressions ([1]) and ([2]), thus the Ra value appears 
to be proportional to m\: 

Ra - 6.64 • \Q'^''m\(j>Alv^f- (H) 

The amount of observed 7-rays that follow the axion ab- 
sorption depends on the number of target nuclei, measure- 
ment time and detector efficiency, while the probability of 
8.4 keV peak observation is determined by the background 
level of the experimental setup. 

3. Experimental setup 

In order to register 8.41 keV 7-rays we used the pla- 
nar Si (Li) detector with the diameter of the sensitive re- 
gion d « 17 mm and 2.5 mm thick. The detector was 
mounted inside a vertical vacuum cryostat « 3 mm below 
the 20 /im thick beryllium window. The Tm203 target 
consisting of 306 mg of ^^^Tm with the diameter of 25 



mm (p=31 mg/cm^) was placed directly on the beryllium 
window. The tantalum collimator of 8-mm radius was sit- 
uated between the beryllium window and detector surface 
to exclude events at the boundary of the sensitive volume. 

The cryostat was enclosed within the concentric copper 
(029 cm) and lead (034 cm) shells, which provided the 
shielding against external radioactivity. The setup was 
located on the ground surface and was assembled of five 
50 X 50 X 12 cm'^ plastic scintillators against the cosmic 
rays and fast neutrons. The rate of 50 ps veto signals was 
600 counts/s, that leads to «3% dead time. The spectrum 
of the Si (Li) signals obtained in the coincidence with veto 
signals allows to check the probability of excitation of the 
8.41 keV level by the nuclear active component and cosmic 
ray muons. 

The data- acquisition system of the Si (Li) detector was 
organized in the following way: the output of the Si(Li) 
detector was fed to the charge-sensitive preamplifier and 
then to two separate shaping amplifiers with different gain 
ratios, which makes it possible to collect spectra from both 
lower (0-;-60 keV) and higher (0-;-500 keV) energy regions. 
This feature allowed us to monitor the background level 
over a wide range of natural radioactivity and, in particu- 
lar, to detect the 7-rays corresponding to the de-exitation 
of the second exited state of ^^^Tm (118 keV). In such a 
way we can evaluate the fast neutron flux which can ex- 
cite 8.41 keV level. The signal from each amplifier was 
received by the individual 4096-channel ADC. Taking into 
account the spectra measured in coincidence with the veto 
system for each amplifier, four 4096-channel spectra were 
recorded. 

The energy scale was defined using standard calibra- 
tion sources of ^^Fe, ^^Co and ^"^^Am. The prompt energy 
resolution of the detector determined by the 14.4 keV 7- 
hne from ^'^Co turned out to be cr=FWHM/2.35=120 eV. 
The high energy resolution and accurate knowledge of the 
energy scale is very important because the energies of the 
characteristic X-rays of thulium are very close to 8.41 keV. 
The most intense L-lines in the case of vacancy on K-shell 
possesses the following energies and intensities: 7.18 keV 
(8.1%, Lai), 8.10 keV (5.2%, L^i) and 8.47 keV (1.6 %, 

L/32)- 

The sensitive volume and the area of Si(Li) detector 
were measured using the X-rays from ^^^Am and ^^Tm 
sources. The self-absorption of 8.41 keV 7-rays by the tar- 
get was found via a detailed M-C simulation. The results 
of simulation were tested with ^"'^Am source placed be- 
hind the Tm203 target. The overall detection efficiency 
for 8.41 keV 7 is estimated to be (1.28±0.06)%. 

4. Results 

The measurements were carried out during 76.5 days 
of live time by brief 2-hour runs in order to monitor the 
time stability of the data- acquisition system. The obtained 
energy spectra of the Si(Li) detector in the range (1 ^ 62) 
keV are shown in Fig[21 While the spectrum measured in 
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Figure 2: The Si(Li)-detector spectra measured during 76.5 days: in 
the anticoincidence (main) and in coincidence (inset) with signal of 
veto system. 



the coincidence with signals of the veto scintillators does 
not contain any peculiarity, one can clearly identify several 
peaks in the spectrum of uncorrelated events. 

Two intense peaks with the energies 13.95 keV and 17.8 
keV are L X-rays of neptunium that is produced in the 

Am a-decay: ^41 Am ^237 rpj^^ ^3 g ^^^y. 

consists of two lines with the energies 13.946 keV (13%, 
Lai) and 13.761 (1.4%, La2), the 17.8 keV peak is more 
complex one formed by L^i_5 lines. The less intense peaks 
with energies of 10.84 keV and 13.2 keV (Lai and L^(i_5) 
of Bi, respectively) are present due to ^lOpb^^iOgj ^_ 
decay in ^ssy series. The well known 59.54 keV Am and 
46.54 keV ^lOpb 7-lines, as weh as 50.7 keV K„i of Tm 
and 57.5 keV K^i of Ta X-rays, were observed in the high 
energy part of the spectrum. The 13.9 keV peak together 
with 46.54 keV and 59.54 keV 7-peaks were used to find 
the final energy scale and the energy resolution a{E) of 
the detector. Since no special electronic stabilization was 
used, the energy resolution determined for the 13.9 keV 
peak(LQi of Np with small admixture of lja2) spread up 
to (7 = 130 eV during the long-time measurements. 

Figl3] shows the detailed energy spectrum within the 
(7.6-^11.4) keV interval, where the "axion peak" was ex- 
pected. There is no pronounced peak at 8.41 keV. In order 
to determine the intensity of the 8.41 keV peak we used 
the maximum likelihood method. The likelihood function 
was determined as a sum of three Gaussians and the linear 



Figure 3: Energy spectrum measured with the Si(Li)-detector in the 
region 7.6-11.4 keV. The fit is shown by sohd hne. 



background assuming that the number of counts in each 
channel had normal distribution. The first Gaussian rep- 
resents the known characteristic Ld X-rays of Bi (Ei = 
10.8 keV, the second Gaussian describes the shape of the 
peak with the energy E2 = 9.2 keV, that we explain by 
Lqi x-rays of Ir (L^i 10.72 keV x-rays of Ir spread 10.8 keV 
peak) and the third Gaussian (Sa) stands for the expected 
8.41 keV axion peak: 



N{E) = a+b-E+ 



1 



271(7 



jexp 



{E.-Ef 



2cr2 



.(12) 



The energy resolution {a) and the position of 8.41 peak 
was fixed during the fitting, while the positions (£'i,i?2) 
and intensities {81,82, 8a) and background coefficients 
(a,&) were independent free parameters. 

The fitting result is given in Fig[3] (solid line). The 
minimum of ~x^/n.d.f. = 121.9/115 corresponds to the 
nonphysical value of the axion peak area 8a = — 14 events. 
The upper limit on the number of events within the peak 
was found via a conventional approach: the dependence of 

on the peak area 8a was calculated for various values of 
8 a while the rest of the parameters were free. Then the ap- 
pearance probability of the given x^{8a) value was found 
and the obtained function P{'x^{8a)) was normalized to 
unity for the S'a > region. The upper limit appeared to 
be equal to 8iim — 31 events for 90% c.l. 

For the given rate of axion absorption Ra the expected 
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Figure 4: The limits on oat, coupling constant obtained by Solax 
[3, Cosme [33, DAMA |3, Tokyo helioscope |i3 and CAST [Til 
experiments. The areas of excluded values are located above the 
corresponding lines 



number of registered 8.41 keV 7-quanta is: 

Sa = e-ryiVi69T™-r-i?A = 3.49-1023.i?^ < 5,™, (13) 

where iViggTm = 1-09 • 10^^ - the number of ^^^Tm nuclei, 
T = 6.61 • 10*^ s - time of measurement, e = 1.28 • 10"^ _ 
detection efhciency and 77 = 3.79-10"^ - internal conversion 
ratio. 

The relation obtained in the experiment (Ra < 
8.8- 10"^"^) limits the region of possible values of the cou- 
pling constant (7^7, 9\n^ 9 an ^'^'^ axion mass ttia- In 
accordance with equations (fTU|) . (fTT|) and on condi- 
tion that (pA/p-y)^ == 1 provided for itia < 2 keV one can 



obtain: 

ffA-y- KffAJV+ffAJvl < 9.2 -10-13 

gA-yiGeY"^) ■ mA{eY) < 1.36 • 10"^ 
ruA < 191 cV, 



(14) 
(15) 
(16) 



aU at 90% c.l. 

The limit ITU) is a model independent one on axion- 
photon and axion-nucleon couplings. The result (jlSp (the 
dimensionless quantity 5^7 w a < 1.36 -10"^^) presented as 
a restriction on the range of possible values of gA-y and niA 
(extracted from gAN) allows us to compare our result with 
the results of experiments using the conversion of axions 
to photon in the laboratory magnetic field or in the field 



Figure 5: The limit on gA^(GeV~^)-mA(eV~^) versus the value of S 
parameter (/3=1, rj=0.5, 2=0.56, solid line ). The doted, dash-doted 
and dashed lines correspond to the values of /3, ij and 2 changed on 
±10%, correspondingly. 



of crystals [|-[2ll. The comparisons are shown in Fig[31 
One can see that our results exclude a region of relatively 
large values of gA-y and rriA ■ For tua ~ 1 keV the limit on 
the gA'y ~ 10"^ GeV"* is only several times lower than 
obtained for low axion masses. 

The limit on the hadronic axion mass {tua <191 eV), 
obtained with the new method, is the strongest one up 
to date among the results of axion detection experiments 
based on the resonant absorption of monochromatic ax- 
ions 24 1- 31 1 or axio-electric effect [i^li^. The obtained 
limit depends on the exact values of the parameters S, 
P, rj and z. As mentioned above and in contrast with 
the 14.4 keV ^^Fe solar axions [3^, the uncertainty of 
the flavor-singlet axial-vector matrix element S does not 
change the obtained constraints significantly: tua <210 
eV (5=0.3) and <180 eV (5=0.7). The changes of 
the nuclear-structure dependent terms from 0.5 to 1.5 
and 77 from -1 to -1-1 lead to the limits niA <(223-;-158) 
eV and ruA < (220-;- 180) eV, correspondingly. The usually 
accepted value of u- and d quark-mass ratio (z= 0.56) can 
vary in the range 0.3-^0.6 [33|, in this case our obtained 
results become tua < (163-^196) eV. The dependencies of 
the limits on gA-yrUA vs S for the various values of /?, rj 
and z are given in Fig[5l 

The sensitivity of the experiment depends on the total 
efficiency of registration defined by the product 77 • e which 



5 



is w 5 • 10~^ in our case. This value can be increased sig- 
nificantly by introducing the Tm target inside the volume 
of a scintillator. 

5. Conclusion 

The search for the resonant absorption of solar axions 
produced by Primakoff conversion was performed. For 
that purpose we used the low-background setup consisting 
of Si(Li)-detector, ^^^Tm target, active and passive shield- 
ing. The obtained upper limit on the values of the coupling 
constant and the axion mass is gA-y{GeV~^) ■ mA{eV) < 
1.36 • 10"*^, which allowed us to set the upper limit on 
hadronic axion mass itia < 191 eV (90% C.L.) (3=0.5). 

The work of E.V.Unzhakov was supported by a grant of 
Saint-Petersburg Government, project no. 2.4/29-04/17. 
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